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We present a new phenomenological model of the dipole scattering amplitude to demonstrate that 
the RHIC data for hadron production in d-Au collisions for all available rapidities are compatible 
with geometric scaling, just like the small- a: inclusive DIS data. A detailed comparison with earlier 
geometric scaling violating models of the dipole scattering amplitude in terms of an anomalous 
dimension 7 is made. In order to establish whether the geometric scaling violations expected from 
small- 2: evolution equations are present in the data a much larger range in transverse momentum 
and rapidity must be probed. Predictions for hadron production in p-Ph and p-p collisions at LHC 
. . . are given. We point out that the fall-off of the transverse momentum distribution at LHC is a 

sensitive probe of the variation of 7 in a region where x is much smaller than at RHIC. In this way, 
, the expectation for the rise of 7 from small- a: evolution can be tested. 
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Cn \ I. INTRODUCTION 

. The observed phenomenon of geometric scaling, i.e. the property that the small- a; DIS cross section depends 

Qh' only on x and through the combination Q^/Q^{x), where Qs{x) is referred to as the saturation scale, still 
requires a satisfactory explanation. It is a property that appears in a natural way in the (asymptotic) solutions 
of nonlinear evolution equations, such as the GLR equation [l], 0] or the BK equation [1, that are expected to 
r*^ become relevant at small x. Nevertheless, the question remains whether the observed DIS data are obtained at 
sufficiently small x values for such evolution equations to be applicable. Here we want to investigate an extension 
of this question to the RHIC d-Au data and the future LHC p-Pb and p-p data. Like the DIS cross section, the 
^ . hadron production cross sections in nucleon-nucleus scattering at high energies have been expressed in terms of 
ps| ' the scattering of a color dipole off small-x partons, which are predominantly gluons 0, [1]. This dipole scattering 
amplitude is the quantity that is expected to display geometric scaling^ and therefore should be a function of 
CO ' Q"^ IQlix)- In nucleon-nucleus collisions the role of Q is played by the transverse momentum qt of the produced 
\ parton that fragments into the observed hadron (or by the inverse of its Fourier conjugate rt). For earlier works 
about d-Au collisions and saturation physics we refer to Refs. 0,3 BUS EI and the review [l^ . 

A successful phenomenological study of experimental DIS data using a model for the dipole cross section was 
. performed by Golec-Biernat and Wiisthoff (GBW) [l^ . They found that the HERA data on the structure function 
F2 at low X (x<0.01) could be described well by a dipole cross section of the form a = aoNGBw{Tt,x), where 
ctq — 23 mb and the scattering amplitude Nqew is given by 
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NGBwirt:X) = l-exp\^~^r^Qlix)j. (1) 

This amplitude depends on x and rt (the transverse size of the dipole) only through the combination rtQ^{x), 
which means it is geometrically scaling. The x-dependence of the saturation scale is given by 



(x) = lGeV {^y^\ (2) 
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^ More precisely, in momentum space times the scattering amplitude is the scaling dimcnsionless quantity. 
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with — 3 X 10 and A ~ 0.3. For nuclear targets contains an additional factor A^^^. 

Model independent analyses of the HERA data show that the low- a; data display geometric scaling for all 

Bin, 

even though the GBW model, (IT])-(I2|jWas found to be inconsistent with newer, more accurate data 
at large and requires modification. In Rcf. such a modification was proposed which includes DGLAP 
evolution, in order to fit the > 20 GeV^ data. In Ref. fvf\ the impact of DGLAP evolution on a geometric 
scaling solution has been numerically studied. An initial condition was constructed so that at — Q'j.{x), with 
Qs{x) as in Eq. the dipole cross section at leading order (and hence as{Q^)xg{x,Q'^)/Q'^) is a constant as 
required for a geometric scaling solution. It was found that under DGLAP evolution to higher values of 
geometric scaling is not violated for A > ANcas/ir in the fixed coupling constant case and only mildly violated for 
all values of A in the running coupling constant case. In the latter case geometric scaling (GS) holds to very good 
approximation in the region log Q^/Q^ <C logQf/A^ (here A denotes Aqcd)- One can conclude that although 
the DGLAP evolution equation for the gluon distribution does not necessarily lead to a GS solution itself, this 
GS property can be preserved to a large extent. The DGLAP induced violations of GS can remain small over a 
wide range of values. 

Similar studies have been performed for the BFKL equation [l^ [l^ . The solution of the BFKL equation with 
an appropriate boundary condition at Qs was found to be geometrically scaling in leading order in the saddle 
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point approximation [lo, [U, 123, llll . Beyond leading order for 1< logQ^/Q^ ^ logQ^/^'^ this solution shows 
approximate scaling. At the scale Qgs = Q'^/A the violations of geometric scaling are considered sizeable. This 
of course assumes that the BFKL equation governs the evolution in this entire region that has been called the 
extended geometric scaling (EGS) region. The region < is referred to as the saturation region. 

The EGS region need not be equal to the region in which GS is observed in experiments, since it is unclear 
that the chosen evolution equation is appropriate in the entire region in the first place. But even if this is the 
case, one could not determine Qgs from the data for a given rapidity. Only if one studies the data as function of 
Q^/Qs fo'' 8. range of rapidities (for which Qgs is not a constant scale) will one be able to establish the extent to 
which GS is violated. 

As mentioned, the DIS data for Q^ > 20 GeV^ prompted the authors of Ref. 16] to propose a modification of 
the GBW model which includes DGLAP evolution. However, in Ref. [23| a model (IIM) has been put forward 
that is a modification of the GBW model and incorporates the violations of geometric scaling expected to arise 
from BFKL evolution in the EGS region. This model leads to a satisfactory fit to DIS data, but without the need 
to include DGLAP evolution at larger Q^. In Rcf. 25] a description of the DIS data is obtained by taking into 
account both BK and DGLAP evolution. The fact that DIS data can be described using different approaches 
suggests that the small-x DIS data do not span a sufficiently large region in Q^ and x to discriminate between the 
different types of evolution. The question is whether or not the RHIC and future LHC data do span a sufficiently 
large region. 

In order to investigate GS violations in the RHIC data, in Refs. [2^ [27l ] a phenomenological model, similar to 
the IIM model, has been put forward (following in part the earlier study of Ref. ^2§\ based on Q). We will refer 
to this model as the DHJ model. It offers a good description of the pt distribution of hadrons produced in d-Au 
collisions at RHIC in the forward region ^, and even in p-p collisions in the very forward rapidity region [29| . 

According to Refs. [1^, [13] the cross section'^ of single-inclusive forward hadron production in high-energy 
nucleon-nucleus collisions is described in terms of the dipole scattering amplitude in the following way. 
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(3) 



A summation over quark flavors q is understood. Here Nf describes a quark scattering off the nucleus, while Na 
applies to a gluon. The parton distribution functions fg/p and the fragmentation functions Di^/g are considered 



^ As it turned out the central-rapidity study of Ref. [27j contained an error in the numerical code. The larger pt data for j/j, = 0, 1 

are in fact not well-described by the DHJ model as will be seen. 
^ To be more precise, Eq. ((Sjl is an expression for the minimum bias invariant yield. 
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at the scale — Pti which we will always take to be larger than 1 GeV^. The momentum fraction of the 
target partons equals X2 = xi exp{—2yh) ■ We find that for pion, kaon, proton and even A production we can to 
good approximation neglect finite mass effects, i.e. we equate the pseudorapidity rj and the rapidity yh and use 
xp = VPt + m'^ / y/s exp{r]) « pt/\/sexp{yh)- Finally, there is an overall iC-factor that effectively accounts for 
NLO corrections. As it is expected that these corrections are more important at small yh, the iiT- factor is allowed 
to be yh dependent. A NLO pQCD analysis of the process pp ^ 7r°X at mid-rapidity for RHIC energies shows 
that such if -factors are relatively constant with pt [30| • 

The dipole scattering amplitude of the DHJ model is given by [1^ H^l : 



Na {quX2) = / d\t e^«* Na {n ,quX2) 



l-exp(^-i(r?Q2(^2))7te.x.)^ 



(4) 



Note that 7 is a function of qt rather than rt ■ This allows one to compute the Fourier transform more easily. The 
corresponding expression Np for quarks is obtained from Na by the replacement (r'^Ql)'' {{Cp /CA)rfQl)^ , 
with Cf/Ca = 4/9. The exponent 7 is usually referred to as the "anomalous dimension", although the connection 
between Na/f a-nd the gluon distribution inside the nucleus cannot always be made. 
The anomalous dimension of the DHJ model is parameterized as 

Ay + + \og{qi/Qi(x2)) 

where y = log l/x2 is minus the rapidity of the target parton. The saturation scale Qs{x2) and the parameter A 
are taken from the GBW model, as given in Eq. Here Qs includes the additional factor A^^^, for which DHJ 
use ^cff — 18.5 in the d-Au case. The parameter d was fitted to the data and set to c? = 1.2. This choice of 7 leads 
to a geometric scaling solution at qt = Qs where 7 = 7s = 0.628 and incorporates to a certain extent the violation 
expected from BFKL evolution for larger qt. The anomalous dimension of DHJ is of the form 7 = 7^ -I- A7, 
where the scaling violations arising from A7 behave as \og{qt/Ql)/y for large y and q| as resulting from the 
analyses of Refs. [2l|, [H, HI] . The question we will address in this paper is whether these violations are really seen 
in the available data. The fact that the DHJ model works well for forward hadron production in d-Au collisions 
does not demonstrate that there are actually violations present as we will show in detail. 
For very large qt, or equivalently for small rt, one can use to good approximation 

iVA(r-t,ft,X2)«i(r2Q2(^2))7te,x.). (6) 

DHJ used the perturbative ^-channel one-gluon exchange result to conclude that 7 — s- 1 as gt — s- cx). However, as 
discussed in Ref. [3l| also if the BFKL equation governs the large-g* region, one can find that 7 — > 1 at large qt. 
The way in which 7 approaches 1 directly determines how fast the cross section will fall off with increasing pt as 
we will discuss in the next section. 

In the DHJ model one retains GS approximately when A7 is small w.r.t. 73. For large, but fixed rapidity 
{y » {d/\Y) this holds numerically up to qt ~ Qgs = Ql/^, when logQ^^/Q^ ~ Xy. In general it is not simply 
the variation of 7 that determines the GS violations. If 7 is chosen to be a function of qt/Ql or r^Q^ only, the 
model is never GS violating no matter how fast it approaches 1 at large qt. We will present such a scaling model 
below and demonstrate that it can describe the RHIC data in both the central and forward rapidity regions. 
Although the new parameterization of 7 is similar in form to that of the DHJ model, it does not have the GS 
violating behavior nor the logarithmic rise expected from the BFKL (and more generally, BK) equation. 

The outline of this paper is as follows. In Section II wc discuss the properties of a new phenomenological model 
of the dipole cross section and the parameters fitted to RHIC data. A comparison to DIS data and the GBW 
model is also made. In Section III we present predictions for hadron and jet production at LHC. In Section IV 
we summarize our main conclusions. 



II. NEW MODEL 



We will now use RHIC data and the cross section expression in Eq. ([3]), as employed by DHJ [H, [23|, to 
constrain the anomalous dimension 7 entering the dipole cross section Especially, we would like to address 



4 



the question whether the RHIC data reaUy require violation of geometric scahng as it was claimed in Refs. |26l.l27| 
and also stated in Ref. [32]. In DIS a scaling behavior of the dipole scattering amplitude N{rQs) maps directly 
into a scaling of the DIS cross section a^.p, which is clearly observable in the data at small x (a; < 0.01). Due 
to the convolution in Eq. ^ the situation is more involved in hadron-hadron collisions where no such scaling in 
terms of the observed kinematic variables {yh and pj) can be expected. Therefore we have to focus on the question 
whether hadron production in d-Au collisions is describable in terms of scaling dipole cross sections Na and Np. 

An anomalous dimension 7 leading to a geometric scaling dipole cross section should depend on w ~ qt/Qs{x2), 
but not separately on qt and rapidity y = logl/a;2. We will take a value 71 = 7(w = 1) of the order of 7^ and 
7^1 for larger w. The parameterization that we adopted reads 

7H=7i + (l-7i) /„ (7) 
(w" — 1) + 

The two free parameters a and b will be fitted to the RHIC data. The parameterizations ([5]) and ([7]) differ not 
only in the scaling behavior, but also the large qt limit of 7 1 is approached much faster in the latter case. 
This will lead to different large momentum slopes of the dipole scattering amplitude ([U and therefore to different 
predictions for the large pt slope using Eq. ([3|). For larger w ~ qt/Qs the exponent can be expanded and the 
dipole scattering amplitude ^ simplifies, cf. Eq. ^ (we will suppress the possible y dependence), 

= _r(-7(u.)) 

for 7 of Eq. (P 



for 7 of Eq. ^ 



For a constant 7 < 1 the amplitude will drop even more slowly than both these models, namely oc / qp'^'^ ■ For 
7 = 1 (the GBW model) one finds on the other hand an unreaHstic exponential fall-off (x exp(— /Q^)/Q^, which 
could be corrected by including a GS violating logarithm as in the MV model [111 . Due to the convolution in Eq. 

with the parton distribution and fragmentation functions, the slope of the pt distribution is not so simple to 
estimate. Empirically we find that the power of the pt distribution is roughly a factor of one to two larger than 
the power of the dipole scattering amplitude. Below we are going to determine this power. We emphasize that the 
fall-off with Pt is not determined by the size of the scaling violations. In order to observe such violations one has 
to study both the yh and pt dependence over a significantly large range. Moreover, the scaling properties of the 
dipole scattering amplitude are not directly visible in the hadron production data, due to the parton distribution 
and fragmentation functions. 



A. Comparison with RHIC data 

In Fig. [T] we show our estimate for dNh/ [dyhd^pt) that follows from the integral in Eq. ^ with our parameter- 
ization for 7(w) (O, which enters the dip ole scattering amplitude All pt distributions of produced hadrons 
measured at RHIC in d-Au collisions |34l . Issl . [36l | are well described. At the saturation scale we have chosen here 
for 7 the same value 7^ = 0.628 as in the DHJ model. We also take A^s — 18.5. We obtain the best fit of the 
data for: 

a = 2.82 and = 168 . (9) 

As mentioned, this LO analysis requires the inclusion of a iiT-factor to account for NLO corrections, which are 
expected to become more relevant towards central rapidity. Following DHJ, the if-factor is allowed to vary with 
but is demanded to be pt independent. The K factors we obtain for yh — 0, 1, 2.2, 3.2, 4 are for our new model 
equal to if = 3.4,2.9,2.0,1.6,0.7 and for the DHJ model K = 4.3,3.3,2.3,1.7,0.7. We have assumed isospin 
invariance to obtain the parton distributions for a deuteron from those for a proton, using the CTEQ5-L0 ones 
[37l |. Furthermore, we use the KKP fragmentation functions of Ref. [38l |. 

From this analysis we can conclude that a GS dipole scattering amplitude is completely compatible with the 
data and therefore the conclusion that GS violations are observed at RHIC cannot be drawn. Of course, a scaling 
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FIG. 1: Transverse momentum distribution of produced hadrons in d-Au collisions as measured at RHIC (black symbols) 
for various rapidities yh- Using the scaling parameterization ^ the data are well described by the expression Q with an 
appropriate if -factor (red/dark curves). The DHJ model © works well only for smaller pt (green/light curves). To make 
the plot clearer, the data and the curves for yh ~ 0,1 and 2.2 are multiplied with arbitrary factors, namely 16, 4 and 2, 
respectively. The STAR data at yh = are from Ref. ^ and 2/^=4 from The BRAHMS results for yh = 1 - 3.2 
can be found in [35l |. 



violating amplitude, i.e. a 7 that depends on w and the rapidity y exphcitly, is not ruled out by the data either. 
What can be concluded further is that the logarithmic rise of 7 resulting from the BFKL evolution incorporated 
in the DHJ model is ruled out in the central region, see Fig. [TJ This may simply indicate that X2 is already so 
large that one is in the DGLAP region. In Fig. [21 the kinematic region where X2 is small is indicated in terms of 
the observables pt and i/h- Where the DHJ model starts to deviate from the data X2 becomes larger than 0.01, 
although Qs is still larger than in DIS at x = 0.01. If one were to exclude the central rapidity RHIC data in the 
model fit, one could also obtain a scaling model with a logarithmically rising, or even constant, 7. 

To indicate how much 7 is constrained by the RHIC data. Fig. [3] shows various 7(u')'s that describe the available 
data equally well. They are all parameterized as in Eq. ([7]) with different a and b values, but require different K 
factors. Clearly, 7 is less well determined close to the saturation scale than in the dilute region. This is because 
the integrand entering the dipole scattering amplitude (jl]) around the saturation scale r = l/Qg is only weakly 
dependent on 7. In addition, the forward data {yh = 3.2 and 4) are essentially sensitive only to 71, since they 
probe the region where w is close to 1. Therefore, the rise of 7 with w is effectively constrained only by the data 
for yh = 0, 1. 

It is important to realize that given a non-scaling "f(w, y) that fits the data for some value of yh one can always 
find a scaling 7(10) that leads to the same pt distribution. This may not be obvious since even if yh is fixed, 
a range of y values is probed in the convolution integral However, the scaling parameter w can always be 
expressed as a function of y and yh, 

'^="~7T= X'i^^viVh]^ = exp[-z/ + y/i]— = w{y;yh,s). (10) 

XF Qs Qs Qs{y) 

Hence, if yh is kept fixed one can express the rapidity y in terms of w and define a scaling "/{w) = 7(w, y — y{w)) 
that leads to the same results as "f{w,y). Clearly, without probing a sufficiently large range of yh scaling violations 
cannot be established. 

As mentioned before, 7 is chosen to be a function of qt rather than vt- For completeness it should be said that 
it is possible to describe the data equally well with a scaling 7 that depends on TtQs- In general, this will be a 
different function than one would obtain by simply replacing qt with 1/rt in 7 which would lead to unphysical 
oscillations in the dipole scattering amplitude and hence in the hadron production cross section. 
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RHIC LHC 

Vs = 200 GeV, d-Au ^ = 8.8 TeV, p-Pb 




FIG. 2: Illustration of the kinematical ranges relevant for RHIC and LHC. The saturation region is set by the line 
qt = Xi/xppt = Qs{x2)- Since the dominant contribution to ([3]) comes from the region xi close to xf = pt/\/^exp{yh), 
we used for this plot X2 — x\ exp{—2yh) ~ pt/v^exp(— y^) and xi « xf- For d-Au we have taken A^s = 18.5 and for 
p-Pb 20. The curves of constant X2 indicate the regions where small- a; physics may become relevant. 



1 I , I , I , I , I , I , [- 




C I I I I I I I I I I I L 

2 4 6 8 10 12 

w = qt/Qs 



FIG. 3: Various fits of 7(u)), which describe the RHIC data equally well. For choices of 7 outside this range the data are 
less well described. 



B. Compatibility with deep-inelastic scattering 

Since the parameterization ^ of the dipole scattering amplitude uses an anomalous dimension 7 7^ 1, the 
resulting amplitude is quite different from the GBW model. Therefore, it is important to check whether our 
anomalous dimension 7(10) is still compatible with the DIS data. For this we use the following expression for the 
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dipole scattering amplitude 

N^{n,Q,x) = l-exp (^_i(r2g2(a;))^(»=VQVQ?W)j , (11) 

where Qs is given by Eq. ^ and for 7 we use our model, Eqs. ([7]) and 

Following the procedure in [l^, we predict the total cross section cr^-j, — ar + crL by folding the resulting dipole 
cross section cr = ctq -^7 with the perturbatively calculable photon wave function, 




where z is the longitudinal momentum fraction of the quark in the dipole. 

In Fig. [4] we show the small-x HERA data fss l . [4 0. l4lj | in a large kinematic range as a function of r = 
/Q^{x). Following Ref. [l^, we scale the HI data [39j| by a factor 1.05, which is consistent with the normalization 
uncertainty. As can be seen, the data for x < 0.01 depend on x and only through the variable r. In Fig. [H 
we compare these data with the original GBW model and the prediction following from our modified 7 obtained 
from a fit to RHIC data. For both models we have neglected effects from finite quark masses in the photon wave 
function, which break geometric scaling. As a result the cross section of the GBW model overshoots the data at 
small T, i.e. at small Q^. Using the modified 7 this effect and therefore the fitted quark mass is smaller since the 
smaller value of 7 in the saturation region suppresses the cross section. Further details of the small-r behavior 
can be found in e.g. [4^ . In addition, we use a somewhat smaller ctq value in order to obtain a better description 
of the data. No parameters of 7 are tuned. The smaller value of (Jq is forced by the region r « 10 . . . 100 where 
7 is not yet close to one but the effective value of rtQs is already large. Given the normalization uncertainty of 
our model, we do not consider the smaller value of do a- problem. Of course, it would be possible to obtain an 
optimized parameterization of 7 by a simultaneous fit to the RHIC and DIS data. But at this stage we conclude 
that the model, which we constructed to describe the RHIC data, can describe the DIS data equally well as the 
GBW model if one adjusts the additional parameter o-q. 




FIG. 4: The 7*p cross section as a function of the scaling variable r measured at HERA [39l, l40ll4ll| . We compare the data 
with the predictions following from the original GBW model, where 7=1, and the modified version fitted to RHIC data, 
where 7 depends on /Q1{x). 



We end this section with a comment on whether or not the factor [Cf/Ca] should have been included in 
of Eq. (jlip . as was done for Np earlier. In order to compare the DHJ model or our new model with the GBW 
model, it would indeed be better to use Eq. as a model for Np and scale {r'jQl)'^ {{CA/CF)rfQl)^ to 
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obtain Na- A fit to RHIC data would then result in a somewhat different 7. Since this is not done by DHJ and 
we are specifically interested in a comparison to the DHJ model, we follow DHJ's approach. This does however 
obscure the comparison to the GBW model somewhat, since that is a model for Np without the factor Cf/Ca- 
Note that for models with 7 7^ 1 this cannot be accounted for by rescaling Qo, because one does not scale Qs 
that enters in 7. In a future combined fit to RHIC and DIS data one would of course like to avoid this slight 
conceptual discrepancy. 



III. LHC PREDICTIONS 



A. Hadron production 



We have seen that where the DHJ model curves deviate from the RHIC data, the a;2-values probed are not very 
small. However, at LHC due to the much higher energies, the region of small X2 extends to a much larger range 
of pt, so that the predictions of the DHJ model and the new model will be different even at small X2- In Fig. [2] 
the region of small X2 is depicted in terms of pt and i/h for p-Pb collisions at LHC. In this section we discuss the 
predictions following from the DHJ model and our new scaling model for the hadron production cross section for 
p-p collisions at ^/s = 14 TeV and for p-Pb collisions at y/s = 8.8 TeV. 

In Fig. [5] we show the predictions for the p-p collisions at y/s = 14 TeV. For smaller pt the predictions of the 
DHJ model and the new model are comparable. This can be expected since this region corresponds to the small-pt 
region at RHIC. For larger rapidities, i.e. yh ~ 7 — 8, the predictions are indistinguishable since the reachable 
momenta qt < ^/sexp{—yh) are so small that the ratios w = qt/Qs are always so close to one that 7 is effectively 
equal to 7^. However, there is quite a large range where the probed values of X2 Pt/^/s exp(— y/i) are small but 
the predictions are clearly different. The slope of the cross section is much larger when described in our model 
as compared with the DHJ model, since 7 rises towards 1 much faster. Hence, a measurement of the slopes at 
moderate rapidities yh at LHC would allow a discrimination between the DHJ model and our model in a region 
where small-x physics may be expected to be applicable. Since a logarithmic rise of 7 is a generic signature of 
BFKL evolution, these measurements offer the possibility of testing whether such small- a; evolution is actually 
relevant at present-day hadron colliders. 

The p-Pb predictions for LHC are very similar. However, due to the smaller energy of ^/s = 8.8 TeV the 
predictions are already comparable for smaller rapidities, i.e. for y^, « 6, cf. Fig. [HI Here the rapidities are given 
for the nucleon-nucleon center of mass frame, which for LHC is not the lab frame in contrast to RHIC. This means 
that in terms of rapidities in the lab frame there is a slight offset of Ay^ = yiab — 2/cm ~ 0.47 to take into account. 

Note that for the whole kinematic range depicted in Figs. [5] and [H] the X2 values are well below 0.01. 



B. Jet production 



Unlike in the case of the DIS cross section, geometric scaling of the dipole amplitude does not lead to scaling of 
the hadron production cross section at RHIC or LHC, because of the convolution of the amplitude with the non- 
scaling parton distributions and fragmentation functions. This effect can be reduced by considering jet production. 
The description of the jet cross section does not involve any fragmentation functions, but reduces to just a sum 
over products of dipole amplitudes and parton distribution functions, 



K{yh) 



dyhtPpt (27r) 



. 9 



fq/pixF,Pt) Np {pt,X2) + fg/p{xF,Pt) Na {pt,X2) 



(13) 



where xp = Pt/\fs exp{yh,) and X2 = xp exp{~2yh) — Ptl\fs exp(— y/j). This means that in the kinematical 
regions where either the gluon contribution or the quark contribution is dominant (in general in a small kinematical 
region), the corresponding distribution function can be divided out, so that one obtains the dipole amplitude 
directly from the data. Of course, if the dipole amplitude is only mildly scaling violating, this kinematical region 
may be too small to observe the violations in this way. At LHC, the gluon contribution to the jet cross section 
is reasonably dominant for transverse momenta pt ^ 15 GeV and hadron rapidities y/i = — 2. In this region, 
the scaling violations of the ratio [p^ dNh/dyh<Ppt)/ fg/p{xF,p1) are in the DHJ model about 30%, while the 
violations for the exactly scaling model ([7]) are, due to quark contributions, still about 10%. We conclude that 
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yh — 0, new model (black), DHJ (green) 
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FIG. 5: Predictions of the transverse momentum distributions of produced hadrons in p-p collisions at the LHC energy of 
^/s — 14 TeV and various rapidities = — 8. The distributions from the scaling model are represented by the black 
lines and those from the DHJ model by the green/light ones. 



pPb hX 



|,V4 ^ = 8.8 TeV 



yii — 0, new model (black), DHJ (green) 

yi, = 2 

yi, = 4 

Vh = 5 

yh = 6 

Vh = 7 
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FIG. 6: Same as Fig. [5] but for p-Pb collisions at the LHC energy of ^/s — 8.8 TeV. We have used AcB — 20. 



it may be difficult to attribute any observed violations directly to Na- A similar conclusion holds for Np in the 
region where quarks dominate (when ^ 0.1, see Fig. [2]) . 

In summary, even for jet production, where there are no complications from the fragmentation functions, it 
may not be possible to establish geometric scaling violations conclusively due to the mixture of quark and gluon 
contributions. The kinematic range at LHC where either quark or gluons dominate is probably too small to reach 
a definite conclusion about scaling violations. 
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IV. CONCLUSIONS 

We have presented a new phenomenological model of the dipole scattering amplitude to demonstrate that the 
RHIC data for hadron production in d-Au collisions for all available rapidities are compatible with geometric 

scaling. Moreover, the model also provides a reasonable description of the small- a; DIS data. On the other hand, 
in a region of yh and pt for which the probed values of x are sufficiently small, the RHIC data are also compatible 
with geometric scaling violating models, such as the DHJ model. The fact that the DHJ model, which incorporates 
scaling violations from BFKL (or more generally BK) evolution to some extent, also describes the forward RHIC 
data suggests that the data simply do not span a sufficiently large region in pt and yh to demonstrate possible 
violations of geometric scaling. Hence, it cannot be concluded that scaling violations of the dipole scattering 
amplitude play a role at RHIC. 

The breakdown of the DHJ model at midrapidity might simply be due to the probed values of x being not 
sufficiently small. The situation is difi^erent at LHC in p-p and p-Pb collisions. For smaller rapidities, but still 
within the region where the small-x description could be applicable, the DHJ model and the new scaling model lead 
to different predictions for the pt fall-off of the cross section. This fall-off is determined by how fast the anomalous 
dimension 7 approaches 1 for large transverse momentum. BFKL evolution typically leads to a logarithmic rise of 
7 — » 1 with transverse momentum and therefore implies a fall-off that is much slower than one finds for the new 
scaling model that is compatible with both the RHIC and the DIS data. Therefore, at LHC in both p-p and p-Pb 
collisions the transverse momentum distribution will probe for the first time at sufficiently small x the rise of the 
anomalous dimension 7, and will thereby provide an important test of the expectations from small-x evolution. 
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